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Abstract 

The Sturmian expansion of the generalized Dirac-Coulomb Green function 
[R. Szmytkowski, J. Phys. B 30 (1997) 825; erratum 30 (1997) 2747] is exploited to 
derive closed-form expressions for electric (<te) and magnetic (<tm) dipole shielding constants 
for the ground state of the relativistic hydrogen-like atom with a point-like and spinless 
nucleus of charge Ze. It is found that cte = Z~ x (as it should be) and 

cr M = -(2Za 2 /27)(4 7l +6 7l 2 - 7 7 i - 12)/[ 7l ( 7l + l)(2 7l - 1)], 



where 7 i = y/1 — (Za) 2 (a is the fine-structure constant). This expression for ayi agrees with 
earlier findings of several other authors, obtained with the use of other analytical techniques, 
and is elementary compared to an alternative one presented recently by Cheng et al. [J. 
Chem. Phys. 130 (2009) 144102], which involves an infinite series of ratios of the Euler's 
gamma functions. 
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1 Introduction 

In the series of papers published by our group over the period of past several years, it has been 
shown that the Sturmian expansion of the generalized (or reduced) Dirac-Coulomb Green function 
(GDCGF), found in Ref. [1], may be used as a convenient tool in perturbation-theory calculations 
of some electromagnetic properties of relativistic one-electron atoms. In particular, closed-form 
expressions, in terms of the generalized hypergeometric function 3.F2 with the unit argument, have 
been derived for the ground-state static dipole magnetizability [2], the polarizability [3] and the 
Stark-induced magnetic anapole moment [4] for the system. 

Recently, Cheng et al. [5] have reported the use of the GDCGF Sturmian expansion technique of 
Ref. [1] for the purpose to find the magnetic dipole shielding constant ctm for the Dirac one-electron 
atom in its ground state. Absolutely no details of calculations have been provided in Ref. [5]; 
only the final expression for ctm has been given therein as a sum of two contributions, one being 
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elementary and the second one having a form of an infinite series of ratios of the Euler's gamma 
functions. A literature search shows that calculations of <tm for the same system were carried 
out before by several research groups and published in Refs. [6-11] (none of those works has been 
referenced in Ref. [5]). An expression for ctm arrived at in Refs. [9-11] (a corresponding formula 
given in Refs. [6-8] contains a misprint, cf. Sec. 3 below) appears to be elementary compared to the 
one in Ref. [5]. This prompts one to inquire whether the GDCGF Sturmian expansion technique is 
practically capable to provide the same simple representation of ctm as the one given in Refs. [9-11] 
(and, after a due correction is made, also in Refs. [6-8]). In the present paper, we answer this 
question affirmatively. 

When preparing this report, we have decided not to focus only on the dipole magnetic shielding 
constant ctm for the Dirac one-electron atom in the ground state, but to present at first in Sec. 2 
details of calculations of the electric dipole shielding constant cr^ f° r the same system. The value of 
the latter quantity is known exactly to be Z^ 1 [13], where Z is the nuclear charge in the units of e. 
This fact makes the evaluation of cte an ideal test of correctness and robustness of any analytical 
technique, including the present one. We believe the material of Sec. 2 is highly instructive, as 
it shows how certain infinite series encountered in calculations based on the GDCGF Sturmian 
expansion technique may be summed to closed elementary forms. The practical knowledge gained 
in that way is then successfully exploited in calculations of ctm reported in Sec. 3. 

A small part of the material of Sec. 3 has been presented in an unpublished comment [12] on 
Ref. [5]. 



2 The electric dipole shielding constant 

Consider a Dirac one-electron atom with an infinitely heavy, point-like and spinless nucleus of 
charge Ze. In the presence of a weak, static, uniform electric field E C xt, an electronic wave 
function of a ground quasi-bound state is, to the first order in the perturbing field, approximated 

by 

tf(r) ~ * (0) (r) + ¥ 1) (r). (2.1) 
Here, ty(°\r) is the ground-state wave function of an isolated atom and is given by 

¥°\r)=a^ 2 (r)+a% 2 *% 2 (r), (2.2) 

where 

K°/ 2 | 2 + |a%| 2 = 1 ( 2 - 3 ) 

(otherwise the coefficients a±\/ 2 are arbitrary), 
with 



= \/ £ r7^T7T (— ) «»<- & /«»>' (2-6) 

y a I (271 + 1) \ a J 

and with r2 KM (n r ), (n r = r/r), being the orthonormal spherical spinors defined as in Ref. [14]. 
The second term on the right-hand side of Eq. (2.1), ^^\r), is the first-order perturbation-theory 
correction to \&(°)(r) given by 

*W(r) = -e£ e xf / d 3 r'G (0) (r,r')r'*(°)(r'), (2.7) 
Jr 3 
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with &°)(r 7 r') being the generalized Dirac-Coulomb Green function associated with the ground- 
state hydrogenic energy level £ ^ — mc 2 ji. 

The symbol j K , appearing in Eqs. (2.5), (2.6), and also in later considerations, is standardly 
defined as 

lK = v/« 2 - (Zar, (2.8) 

where a is the Sommerfeld's fine-structure constant. Moreover, as usual, ao denotes the Bohr 
radius. 

The electric field produced at the point r by the atomic electron, the latter being in the state 
characterized by the wave function \I/(r), is given by 



E 

where 



(r) = -if dV Pc (r')^4, (2.9) 
4ire J W 3 \r - r'\ 3 



p (r) = -e&(r)^(r) (2.10) 



is the smeared electronic charge density distribution. Hence, at the point r — 0, where the nucleus 
is located, the field is 

J3(0) = -J- / dVr'^. (2.11) 
47reo J R 3 r' 3 

Now, if \E"(r) is approximated as in Eq. (2.1), we have 

p c (r)^p(°)(r) + pW(r), (2.12) 

where 

pi°\r) = -e*W(r)*M(r) (2.13) 

and 

pi^ir) = -2eRe[* (0)t (r)*( 1 )(r)]. (2.14) 
Consequently, E(0) is approximately given by 

E(0)~ £ (0) (0)+£ (1) (0), (2.15) 

with 

B ..) W ^f J y/""M?" l W p.,,, 

4ire J R3 r' 3 

and 

g (.) (0) _ * Re f rfV^^W" "<--). (2 . 17) 

It may be easily shown that, because of parity reasons, the field component at r = due to the 
unperturbed electronic distribution, E^°\o), vanishes. 

In turn, using the expression (2.7) for ^^\r), we find that the first-order field correction 
^^'(O) may be written in the form 

E {1 \o) = -S E -E cxU (2.18) 
where Se is the electric dipole shielding tensor given by 
2e 2 



Sr = 



47TC0 

To evaluate S E , we rewrite Eq. (2.19) in the form 
2e 2 



Re/ d 3 r / d 3 r'*(°)t(r)^G< »(r,r')r'1'(»)(r'). (2.19) 

Jr3 J R 3 r 3 



Se = 



47T£o 



Re V e n e* n , f d 3 r f d 3 r'^{r)r- 2 e* n -n r G^(r,r')r'e n >-n' r ^°\r'), (2.20) 



3 



where e n , with n — 0, ±1, are the unit vectors of the cyclic basis, related to the Cartesian unit 
vectors through 

e Q = e z , e±i = ^-j=(e x ± \e v ). (2.21) 

In the next step, we substitute (r), as given by Eqs. (2.2) and (2.4), and the multipolc expansion 
of the generalized Dirac-Coulomb Green function, which is 

, oo \K\-\/2 

G (0) (r,r') = ^ £ E ^7 

A-=-oo M=-\K\+l/2 
(KjiO) 

x / .9^ (++) (^r')^/fMK)f7^ M «) -ig ( ^ (+ _ ) (r,r')n K M(n r )nl KM (n' r ) \ 
\ ijq_ +) ('','-')J!-KMK)^MK) ^!(— ){r,r')n-KM(n r )rt_ KM (n' r ) /' 

(2.22) 

into Eq. (2.20). After the angular integrals appearing in the resulting series representation of Se 
are evaluated with the aid of the following identities obeyed by the spherical spinors [14]: 



2 M J(k+\) 2 -H 2 
e • n r VL Kll (n r ) = - ^ 2 - - Sl- Kli (n r ) H + to K +i,n(n r ) 



2 - M 2 



+ V |2w _ ^ r! K -i^(n r ), (2.23) 



Jk 2 - (/i± ±) 2 
e±i • n r O re/J (n r ) = ± V2 ^ 2 _ ^ Q_ KijU ±i(n r ) 

+ V2(2. + l) - n ^W^) 

- V2(2.-l) < 2 - 2 ") 

we arrive at the conclusion that for arbitrary (up to the normalization constraint (2.3)) values of 
the coefficients a^j/ 2 the electric shielding tensor for the hydrogenic ground state is a multiple of 
the unit dyad, 

S E = <7 E I, (2.25) 

and that the only contributions to the shielding constant cte come from the terms with K = +1 
and K = —2. One has 

ce = oe,i + <j e,-2, (2.26) 

where 

= 2 -J™ArJ™dr> ( P(0) (r) Q«>)( r ) ) r 2 G< '(r,r')r' ( ) (2.27) 



and 



with 



<7 E; _ 2 = - jf dr jf dr' ( P(0)(r) Q(°)(r) )r- 2 G^(r,r>' £ (0) £,j J , (2.28) 



G^(rr')-( 9 {0) K A++) {ry) g<® *,(+-) (r,r') \ 



(2.29) 
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being the radial generalized Dirac-Coulomb Green function associated with the hydrogenic ground- 
state energy level. 

To evaluate the integrals in Eqs. (2.27) and (2.28), we shall exploit the following separable 
expansion of G^(r,r'): 



GP(r,r')= £ 



UK ) ( MM&M T%(r>) ) (K * -1), (2.30) 



« (0) - 1 



n=-oo r'nK 

which is a corollary from the theory of the Dirac-Coulomb Sturmian functions [1]. Here 



I (l + 7i)(M+2 7 x-)|n|! 
2ZN nK (N nK - K)T{\n\ + 2 lK ) 



L 



(27k) 

H-i 



2Zr 
a 



K-N nK (2Zr\ 
\n\ + 2 lK M 



- 



(2.31) 



and 



2$(r) 



/ (l- 7 i)(|n| + 2 7 K)|n|! 
2ZN nK (N nK - K)T{\n\ + 2 lK ) 



-Zr/ao 



{2lK) (2Zr\ 



K_NnK L {2- <K ) (*Zr\ 
\n\ + 2 lK M \ a Q J 



(2.32) 

(with L^\p) denoting the generalized Laguerre polynomials [15]; we define L^l(p) = 0) are 
the radial Dirac-Coulomb Sturmian functions associated with the hydrogenic ground-state energy 
level, and 

(0) _ M + IK + N nK 



71 + 1 



(2.33) 



where 



N nK = ±^{\n\+ lK ) 2 + {Zaf = ±v/|n| 2 + 2|n| 7 K + if 2 



(2.34) 

is the 'apparent principal quantum number' (notice that it may assume positive as well as negative 
values!). The following sign convention applies to the definition (2.34): the plus sign should be 
chosen for n > and the minus one for n < 0; for n = one chooses the plus sign if K < and 
the minus sign if K > 0. 

At first, we attack the double integral in the expression for <te,i- Inserting Eqs. (2.5), (2.6) and 
(2.30)-(2.33) into the right-hand side of Eq. (2.27), taking the two resulting separated integrals 
with the aid of the known formula [16, Eq. (7.414.11)] 



Jo 



, A- PrWf;)1 _ r(/3 + l)r(rc + tt-/3) _ rC8 + l)rQ8-a + l) 
app L nKP) - nlT{a _ p) — nW(J3-a-n+l) 



and employing Eq. (2.34), we obtain (7 E ,i in the form of the following finite sum: 

2 



[Re(/3) > -1] 
(2.35) 



CE,1 



9Z 



(7 1 + l) 2 (2 7l + l)r(2 7l -l) 



n=-2 



(N nl - l)(jnj + jV„i + l)(|n| - 7l JV nl -7i + l) 
JV„i(2 - |n|)!r(|n| + 2 7l + l)(|n| + N nl - 1) ' 



(2.36) 



Using again Eq. (2.34), after some algebra, the right-hand side of Eq. (2.36) simplifies considerably, 
yielding 

271 + 1 (9 , 7 x 

^ ~ ~9Z~- (2 ' 37) 
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As one of the two separated integrals evaluated above converges at its lower integration limit 
provided 71 > 1/2, we have the following constraint on the nuclear charge: 

Z < c^ 1 ^- - H8.67. (2.38) 

Evaluation of the double integral on the right-hand side of Eq. (2.28) appears to be much more 
cumbersome. Proceeding initially as in the case discussed above, employing Eq. (2.34) and the 
trivial but useful identity 

72 - 7i + 3, (2.39) 

with much labor we obtain 

1 (7i + l)r(72+7i-l)r(72+7i+2) 



CE.-2 



2Z r( 72 - 7l - i)r( 72 - 71 + 2)r(2 7l + 1) 

r(M + 72 - 71 + i)r(M + 72-71-2) JV n ,_ 2 + 2 



00 

X 

n— — 00 



E ln|!r(|n| + 2 72 + l) N n> _ 2 

(H + 72 + 71 + 1 - 7i^»,-2)[(7i - l)(27i + 5) - (2 7 i - 1KH + 72 + N n - 2 )] 

\n\ + 72 - 71 - 1 + Nn,-2 

(2.40) 

Collecting together those terms in the above series which correspond to the same absolute value 
of n, again with the use of Eqs. (2.34) and (2.39), we arrive at 

4 (7i + l)r(72+7i-l)r(72+7i+2) 

CE,-2 — — 



9Z r( 72 - 7i - i)r( 72 - 7i + 2)r(2 7 i + 1) 

r(n + 72 - 7i - 2)r(n + 72 - 71) 



E 



x [3(71 - l)(n + 72 - 7i - 2)(n + 72 - 71) + ( 2 7i 2 + 271 - 3)(n + 72 - 71) - 3(7i - 2)] . 

(2.41) 

The series in Eq. (2.41) docs not terminate. Nevertheless, we shall show that oe,-2 may be 
expressed in the form which is as elementary as the one in the case of <te,i- To this end, we invoke 
the relationship 

^ r(n + Ql )r(n + fl2) z n r(ai)r(a 2 ) „ / 01,02 \ / M <- n (oao\ 

where 2-F1 is the hypergeometric function. Transforming Eq. (2.41) with the aid of the above 
formula results in 

4 (7i + l)r(72+7i-l)r(72+7i+2) 

CE,-2 



9Z (72 - 71 - 2)( 7 2 - 71) (72 - 7i + l)r(2 7 i + l)r(2 72 + 1) 

72 -71 - 1,72-71 + 1 



3(7i - 1)(72 - 7i - 2) (72 - 7i) 2F1 
+ (27 1 2 + 2 7l -3)(72-7i)2^i 



2 7 2 + l 

72-7i-2,72-7i + l 
272 + I 



o\ p I 72-71-2,72-71 , 
3(71-2)2^) 2 ^ 2 + 1 ,1 



(2.43) 



Now, we owe to Gauss the following identity [15]: 
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Applying it to the three 2 Fi functions appearing in Eq. (2.43), we find 



and 



2-F1 

2 -Pi 



72 - 7i - 1 i 72 - 7i + 1 
272 + 1 

72-71^2,72-71 + 1 
272 + 1 

72 - 7i - 2,72 - 7i 



1 = 



272 + I 



f 



r(2 7l + i)r(2 72 + f ) 
r( 72 + 7i) r (72 + 71 + 2) ' 

r(2 7l + 2)r(2 72 + 1) 

T(72 +7i) r (72 +7i+3) 

r(2 7l + 3)r(2 72 + i) 



(2.45) 
(2.46) 

(2.47) 



r( 72 + 71 + !) r (72 + 71 + 3) 

Plugging Eqs. (2.45)-(2.47) into Eq. (2.43), after some further simplifications based on the use of 
the relation (2.39), we arrive at the afore-announced elementary representation of the K = — 2 
component of the shielding constant: 



CE,-2 = 



2(7i - 4) 
9Z ' 



(2.48) 



Insertion of (Te,i and oe,-2 as given in Eqs. (2.37) and (2.48), respectively, into Eq. (2.26) leads 
us to the conclusion that the electric dipole shielding constant for the Dirac hydrogenic atom in 
its ground state is 

a E = Z-\ (2.49) 



as it should be. 

In the nonrelativistic limit 



so that one has 



7l ~ 1 - -(Zaf, 



and 



ce.i 



CE,-2 



3Z 
2 

iZ 



l-^Zaf 



l + -(Zaf 



(2.50) 
(2.51) 
(2.52) 



3 The magnetic dipole shielding constant 



We proceed to the evaluation of the magnetic dipole shielding constant. The model of the hydrogen- 
like atom we shall adopt for this purpose is the same as in the preceding section. 

In a weak, constant, uniform magnetic field B C xt, the atomic ground energy level splits 
into two, their energies being given, to the first order in -B ex t, by 



E^E^+E^ (M = ±i), 



with 



(3.1) 



(3.2) 



where [Ib is the Bohr magneton. The corresponding wave functions, to the same approximation 
order, are 

^(r)^*W(r) + *«(r) (p = ±|), (3.3) 

with <4 0) (r) given by Eq. (2.4) (the space quantization axis being chosen along the external 
magnetic field direction) and with 



¥V(r) = - l -ecB C xf I d 3 r'GW(r,r')(r' xa)^(r r ). 
2 Jr3 



(3.4) 
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Here, a is the Dirac 4x4 vector matrix given standardly by 



«-(;;). («) 

where <x is the vector composed of the Pauli matrices. 

According to the Dirac theory, there is an electric current with density 

J„( r ) = -ec*j,(r)atf„(r) (3.6) 

associated with the electron being in the state \1/ M (r). Using the Biot-Savart law, we find that at 
the point r the magnetic field due to the current distribution j M (r') is 



In the particular case when the observation point is located at the nucleus, i.e., at r = 0, the above 
expression simplifies to 

,3 / r ' X •?'>') 



(0) = £ / dV r X y° . (3.8) 

If, as we have assumed above, the electronic wave function Vl/ At (r") is known to the first order in the 
perturbing field B cxt , the current j Ar) may be approximated as 

i.M^fW + ^W, (3-9) 

with 

fi\r) = -ec*f\r)cc^\r) (3.10) 

and 

j«(r) = -2ecRe[*(°)t(r)a*«(r)]. (3.11) 
Consequently, for the magnetic field at the nucleus location we have 

B,(0),B(»»(0) + BW(0), (3.12) 

with 

*<«»(0) = / dV jgV>^^ (3 . 13) 

and 

4tt 7 R 3 r' 3 

If Eq. (3.13) is rewritten in the form 

Bf (0) _ -eo£ £ </ E /V «g"W«--W'"')»g'C J ), (1,5) 
using Eqs. (2.4)-(2.6) and the spherical spinor identities [14] 



^Ti fl -«(»K) 1 ' ' 

.V / («-|) 2 -M : 



e • (n r x <x)0 KA1 (n r ) = i ^ 2 - i »_ KAt (n r ) +i i 2k+ i| ^K+i,M( n r) 



|2k- II 



2 

fi«_i )At (n r ), (3.16) 
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/ K 2_ (M± 1)2 

e±i • (n r x er) K([1 (n r ) = =fi2V2k ^ 2 - ^ to- K ,p±i(n r ) 

+ i v^TT) n ^ 1 ^) 

v/(«TM - |)(«TM - |) 
+ 1 V2(2.-l) ( 3 - 17 ) 



with no difficulty one finds that 



where 



B(» = ^M ^^.!) ^, (3-18) 



&o = ^ (3.19) 



47r a; 



is the atomic unit of the magnetic field induction, while n ext is the unit vector along £? e xt- The 
radial integral 



poo 

/ drr- 2 p(°)(r)Q(°)(r), 
Jo 



encountered in the course of evaluation of B^(0), converges at its lower limit provided 71 > 1/2, 
which constrains the nuclear charge as in Eq. (2.38). 

To evaluate the first order approximation to the magnetic field induced at the nucleus, we plug 
Eq. (3.4) into Eq. (3.14). This gives 

B( J 1) (0) = -S M ^-Bcxt, (3.20) 

where 

£ MAJ = --^Re/ d 3 r / dV*( >t(r)^G( )(r,r')r'x«*(°)(r') (3.21) 
47re J R 3 J r3 r J 

is the magnetic dipole shielding tensor. Rewriting Eq. (3.21) in the form 

S M , M = -^-Re V e n e* n , [ d 3 r / d V $> { °^ (r)r~ 2 ■ (n r x a)G (0) (r, r')r'e n / • « x a)*™ (r), 

^0 n y^_ x j ua v 

(3.22) 

inserting the multipole representation (2.22) of the generalized Dirac-Coulomb Green function 
into Eq. (3.22) and using the relations (3.16) and (3.17) to carry out integrations over the angles 
of the vectors r and r', after some labor we discover that, like its electric counterpart, the tensor 
is a multiple of the unit dyad: 

^m,/j, = ct m'- (3.23) 

As the notation used indicates, the factor cxm, being the magnetic shielding constant for the system 
under study, is independent of the angular momentum projection quantum number 11; it is given 
by the sum 

cm = &m-i + cm, 2, (3.24) 



where 

cm, 

and 



1 = dr dr' ( Q(0) (r) P«»( r ) ) r" 2 ^, r>' ( g^j ) (3-25) 



cm, 2 = -\j~drj~dr' ( Q(0)(r) p(°)(r) ) r" 2 ^ (r,r>' ( g^j ) . (3.26) 
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The Sturmian expansion of G { £\r, r') given in Eq. (2.30) is inapplicable in the case of K = — 1 
and thus it cannot be used for the purpose of evaluation of the double integral in Eq. (3.25). The 
valid Sturmian representation of G^(r, r') which we employ for that purpose is [1] 

= £ V-r f $> t] ) ( ) 

n=-oc Mn,-1 1 \ n,-lV/ / 
(njiO) 

+ ("« - i) ( ) ( C(r-) !•»!,(/, ) 

) ( C>M if-M ) 



+ 



Swl *<V)). (3-27) 



lfi(r) 



where 



/W(r) = ( 71 - |) l( r) +71 + Za) l( r), (3.28) 

j(o) (r) = /(0) (r) + 5 W i(r) 

= (7i + i)^°l 1 M+7i( 1 ^f + ^)ri°2 1 (r) (3.29) 



and 



^ (0) (r) = 7i (^f - S$Ur) - (71 - J) ^(r). (3.30) 
Exploiting Eqs. (2.5), (2.6) and (2.31)-(2.35), we find that 



- - ^-(7i + l)(27i + l)r(2 7 i - 1) 



\p f ]n N n ,-i + l \n\ + l 

„_ 2 Ar "." 1 (2-H)!r(|n|+2 7 i + l)(|n|+JV n ,_ 1 -l) 

(n^O) 

2Za 2 2 7l + l 2Za 2 2 7l + l 2Za 2 2 7l + 1 
9 7i 9 2 7 i - 1 9 7i ' [ ' j 

terms on the right-hand side being ordered in correspondence with Eq. (3.27). The first of them 
is readily found to be 8Za 2 /[9 7 i(2 7 i — 1)], so that <7m,-i is 

2Za 2 2 7l 2 + 7l -4 
f ™'- 1 -"^" 7i(27i-l)- (3 ' 32) 

To calculate <tm,2 from Eq. (3.26), we use G 2 °^(r, r') in the form (2.30). Applying Eqs. (2.5), 
(2.6), (2.31)-(2.35) and (2.39), we arrive at 

2Za 2 T( 72 + 7l - l)r( 72 + 7i + 2) 

cm, 2 — 



9 r( 72 - 7 i - i)r( 72 - 7l + 2)r(2 7l + 1) 

r(|n| + 72 - 7i - 2)r(|n| +72-71+2) N n2 - 2 



oo 

X 



E 



M!r(M + 2 72 + i) N n2 

3|n|+3 72 + 7 i + l + 3^„2 (3 33) 

M + 72 - 71 - 1 + N n2 
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The simplification is achieved after one collects terms with the same absolute value of n. This 
gives 

2 Za 2 r( 72 +7i-l)r(72+7i+2) T(n + 72 - 7l - l)T(n + 72 - 7l + 2) 

9 r( 72 - 71 - i)r( 72 - 7l + 2)r(2 7l + 1) f- n 



cm, 2 



n=0 



n!r(n + 2 72 + l)(n + 72 - 7l ) 

(3.34) 

which is the final expression for cm, 2 provided in Ref. [5]. We shall show that, similarly to the 
series in Eq. (2.41), the one in Eq. (3.34) may be also summed to a closed elementary form. To 
this end, we rewrite the latter equation as 

2Za 2 r( 7 2+ 7 i-l)r( 72 + 7 i + 2) 



cm, 2 



E 

OG 

+£ 



r( 72 - 71 - i)r(72 - 7i + 2)r(2 7l + 1) 

T(n + 72 - 7i - l)r(n + 72 - 71) 
nW(n + 2 72 + 1) 

T(n + 72 - 71 - i)r(n + 72-71 + I) 



n\T(n + 2 72 + 1) 
and further, after use is made of Eq. (2.42), as 

2Za 2 r( 72 +71 - l)r( 72 +71 +2) 



(3.35) 



cm, 2 — 



9 (72 - 71X72 - 71 + i)r(2 7 i + i)r(2 72 + 1) 

72 - 71 - 1)72 - 7i 



2*1 



272 



1 



; 1 + (72 - 71) 2-F1 



72-71-1)72-71 + 1 
272 + I 



According to Eq. (2.44), the first of the two 2-F1 functions appearing above is 

r(2 7 i + 2)r(2 72 + 1) 



2F1 



72 - 71 - 1)72 - 71 

272 + 1 



1 = 



(3.36) 



(3.37) 



r(72 + 71 + i)r(72 + 71 + 2) ' 

while the second one has been already encountered in the course of evaluation of <xe,-2 and is given 
by Eq. (2.45). Hence, once Eqs. (2.45) and (3.37) are inserted into Eq. (3.36), we eventually arrive 
at the following simple representation of cm, 2: 

2Za 2 71 + 2 

CTM ' 2 = 27 71 + 1- ( } 

From the expressions for the components cm,-i and % 2 of <7m, displayed in Eqs. (3.32) and 

(3.38) , we deduce that the total magnetic shielding constant for the hydrogen-like atom in its 
ground state is 

2Zc? 4 7l 3 + 6 7l 2 - 7 7 i - 12 
aM = - 27 71(71 + 1X271-1) ' (3 ' 39) 
in agreement with the findings of Moore [9], Pyper and Zhang [10] and also of Ivanov et al. [11] 
(after one takes into account that the latter authors define cm with the opposite sign). The 
expression for ctm given by Zapryagaev et al. in Refs. [6-8] becomes equivalent to the one in Eq. 

(3.39) after in the denominator of the former one replaces A2 by Ai, the latter quantity being 
identical with our 71 . 

In the nonrclativistic limit one has 

a M ,-i ^ [1 + 4(Za) 2 ] (3.40) 



and 



and also 



cm, 2 



cm 



j Za 2 
9 

Za 2 



l+ljW 



(3.41) 
(3.42) 



11 



4 Conclusions 



In this paper, we have proved that the analytical technique based on the Sturmian expansion of the 
generalized Dirac-Coulomb Green function enables one to arrive at the same simple closed-form 
representations of the electric (cte) and magnetic (ctm) dipole shielding constants for the Dirac 
one-electron atom in its ground state as have been given previously by several other authors. We 
find this encouraging in the perspective of our planned, mathematically much more challenging, 
application of the same calculational technique to the evaluation of a M for an arbitrary state of 
that atom. Accomplishment of this task would be undoubtly worthy of effort since although results 
for some particular excited states are known (cf., e.g., Ref. [9-11]), in the most general case, to the 
best of our knowledge, calculations of ctm have never been carried out. 
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